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The simulation of quantum correlations with alternative nonlocal resources, such as classical communication, 
gives a natural way to quantify their nonlocality. While multipartite nonlocal correlations appear to be useful 
resources, very little is known on how to simulate multipartite quantum correlations. We present the first known 
protocol that reproduces 3-partite GHZ correlations with bounded communication: 3 bits in total turn out to be 
sufficient to simulate all equatorial Von Neumann measurements on the 3-partite GHZ state. 
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When measurements are performed on several quantum 
systems in an entangled state, the statistics of the results may 
contain correlations that can't be simulated by shared local 
variables. Such correlations are called nonlocal. They can be 
identified by their capacity to violate some inequality: these 
so-called Bell inequalities are indeed satisfied by all correla- 
tions that can be explained by shared local variables (Q]]. 

The observation that quantum theory predicts nonlocal cor- 
relations is not new; it goes all the way back to the famous 
EPR argument [2J. Many experimental confirmations have 
been demonstrated all over the world during the last two 
decades of last century Q]. During the first ten years of this 
century, the interest for nonlocal correlations has shifted from 
mere skepticism and incredulity to more constructive ques- 
tions. First, physicists raised the question of the power of 
nonlocal correlations for information processing; the main 
examples being "device-independent" quantum key distribu- 
tion 14|-|6|] and random number generation H. Second, the- 
orists realized that quantum correlations, although possibly 
nonlocal, are never maximally nonlocal, hence the question 
"why is quantum theory not more nonlocal?" (Ht; note the 
great advances since the original question "Why is quantum 
theory not local?". 

Thirdly, and this is the topic of this letter, physicists and 
computer scientists tried to quantify nonlocality; that is, to 
treat nonlocality as a physical quantity. Indeed, the violation 
of a Bell inequality only proves that the correlations are not 
local, but doesn't tell us anything about how far from local 
they are, i.e. how much nonlocality they contain. Intuitively, a 
larger violation should signal more nonlocality. But this naive 
approach is insufficient as some correlations may violate dif- 
ferent Bell inequalities by different amounts. A quite natu- 
ral measure of nonlocality is the number of classical bits that 
need to be communicated from one party to another in order to 
simulate the correlation. For local correlations, no communi- 
cation is needed, as shared local variables suffice; hence local 
correlations have a "communication measure" equal to zero, 
as it should be. Let us stress that the idea is not to imagine 
that nature does use communication to produce nonlocal cor- 
relations, it is only to quantify the amount of nonlocality by 
the quantity of communication required to simulate the corre- 
lation. 



That such a measure of nonlocality is natural is testified 
by the fact that it has been introduced by 3 independent pa- 
pers |0- ill- More precisely, in this letter we adopt as a mea- 
sure the number of bits communicated between all partners in 
the worst case J9l LL0|] . An alternative could be to count the 



number of bits sent on average 111 ll 1 1 2 ] . 



For the case of 2 qubits in a maximally entangled state, 
Toner and Bacon H 1 3fl proved that one single bit of communi- 
cation suffices (if one restricts the analyses to projective Von 
Neumann measurements, as we do in this letter). Hence the 
nonlocality of two spin-^ particles in the singlet state is 1 bit. 
For the general case of 2 qubits in a partially entangled state 
it is known that 2 bits of communication are enough 11311 . 
though it is still unproven that one bit isn't sufficient. At 
first, one may think that the nonlocality of a partially entan- 
gled state shouldn't be larger than that of maximally entangled 
states, but this is not so clear once one realized the difficulty 
of simulating at the same time the nonlocal correlation and the 
nontrivial marginal probabilities IU4[|l5ll . 

GHZ correlations. In this letter we consider 3-qubit GHZ 
quantum correlations and present the first known protocol to 
simulate such nonlocal correlations with bounded communi- 
cation. This problem is the straightforward next step after the 
2-qubit case; it attracted the attention of most of the special- 
ists. After years of unsuccessful efforts, the feeling started 
to spread that it might be impossible with finite communi- 
cation flal . Some hope, however, appeared when Bancal et 
al. fl7ll presented a protocol with unbounded, but finite aver- 
age communication. Moreover, a team recently presented a 
nonconstmctive existence proof of a protocol with 6 bits of 
communication ifTill : the proof turned out to be flawed, but 
the impulse was given! 

More precisely, our goal is to simulate, with classical com- 
munication, the quantum correlations obtained by performing 
equatorial Von Neumann measurements on a 3-partite GHZ 
state. Namely: 3 parties, Alice, Bob and Charlie, each receive 
an input angle 4>a,4>b an d 0c £ [0, 2tv] (corresponding to a 
measurement setting on the equator of the Bloch sphere S 2 ), 
and they must output binary outcomes a, (3, 7 e {+1,-1}, 
such that the expectation values satisfy 

(afti) = cos(</>a + §B + 0c), (1) 
while all single- and bi-partite marginals vanish. Note that 
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although the choice of equatorial measurements is restrictive, 
these are enough to come up with the "GHZ paradox" 1119 1. 
We will show that our problem can be solved with finite com- 
munication. For that, we first introduce a protocol that pro- 
vides "stronger" correlations, before showing how to ade- 
quately transform these and obtain the desired cosine corre- 
lations. 

Simulation with classical communication. Consider the 
following protocol, that uses 3 bits of communication: 2 from 
Bob to Alice, and 1 from Charlie to Alice. The sign function 
below is defined as sign(:c) = +1 if x > 0, sign(x) = —1 if 
x < 0. 

Protocol 1. Let Alice and Bob share two random vectors \\ 
and A2, uniformly distributed on the sphere S 2 , together with 
a random bit £ £ {0, 1}; let Alice and Charlie share a random 
variable ip c , uniformly distributed on [0, 2n]. 

After reception of their measurement settings 4>a, 4>b and 
(j>C, the three parties proceed as follows: 

0. Bob defines b to be the equatorial vector with azimuthal 
angle ^ — 20^; he calculates a$ = sign(6 • Ai)sign(6 • 

A2), and sends the bit To = 1 ~ CT ° to Alice. Alice and 
Bob can then both determine the azimuthal angle ipo S 
[0,27r] of Ao = Ai + ( — 1) T °A2," they calculate iff, = 

£tt g [0, 24 
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1. Alice, Bob and Charlie define (j) a = 4>a~ <fb~ fc, 4>B = 
4>B + fb an d 4>c = 4>C + fc, respectively. 

2. Bob calculates Ub = sign(sin20s), and sends t\> = 
1 ~ CTb to Alice; he outputs (3 = sign(sin 

Similarly, Charlie calculates a c = sign(sin 20c), and 
sends t c = 1 ~ 1Tc to Alice; he outputs 7 = sign(sin 0cO- 



3. Alice outputs a = sign(sin(— (f>A — T\pk — t c ^)). 

Before analyzing the correlation given by Protocol 1, let 
us give an intuitive understanding of it. Forget for now the 
rather technical step ll20ll . and note that after step 1, one has 
0A + 0s + 0c = 4>A + 4>b + 4>c ; the first two steps will ensure 
that the final tripartite correlation depends on the sum <fi = 
4>A + 4>b + 4>c only, and that all marginals vanish. Assume 
now that 4>b,4>c £ [0, 7r] (and hence (3 = 7 = +1); if this is 
not the case, Bob and Charlie can locally subtract tt to ipB or 
<f>c and flip their output, so that the correlation E(4>) = (a(3j) 
is unchanged - this is precisely why we ask them to output 
/3 = sign(sin0B) and 7 = sign (sin )■ In step 2, Bob and 
Charlie tell Alice in which quadrant ([0, or [5,7r]) their 
angles <j>B and <j>c are. From this information, Alice knows 
in which half-circle <ps + <t>c is (more precisely, she knows 
sign(sin(0s + 4>c)) ° r sign(cos(0s + ^c))> depending on 
whether t/, = t c or t/, ^ r c ); if —cf>A is in the same half-circle, 
she wants to obtain a good correlation with /3j = +1 (if by 
chance — <\>a = 4>b + <frc, she wants a perfect correlation!), 
and will thus output a = +1; otherwise, she will output a = 
— 1; this corresponds precisely to step 3. 




FIG. 1: Correlation Ei(4>) — Ei(4>a + 4>b + 4>c) obtained from 
Protocol 1, compared to the desired correlation cos(<j!>). 



As shown in Appendix A, Protocol 1 gives vanishing 
marginals, and the following 3-partite correlation E\ (<fi) := 
(«/?7): 

■30 1 1 

M4>) = -,Y,7KT7^ . cob ((2n+lfr) 



tt 2 ^ (2n+l) 2 4-(2n+l) 2 

n >0 



1 



2(f) — sin 20 



for cf> G [0, 7r] 



(2) 



which, as akeady mentioned, only depends on the sum = 

4>A + 4>B + 4>C- 

Ei{4>) is shown on Figure Q] One can notice that it is 
"stronger" than the desired cos <\> correlation, in the sense that 
l-Ei (<f>) I ^ I cos 01 f° r a U 4>- Intuitively, one should be able to 
add some noise and weaken the correlation. However, weak- 
ening any given stronger correlation so as to obtain the desired 
cosine is not trivial, since this weakening must depend on <j> 
and should in particular not weaken the extreme correlations 
for <j> = and tt. In fact, it seems that correlations must in 
general have quite specific properties for them to possibly be 
transformed to the desired cosine. 

In order to do so, and starting from a 27r-periodic correla- 
tion function such that E(0) = —E(ir) = 1, one can for in- 
stance try to mix correlations of the form E((2m-\- 1)0), with 
to G Z , as this will preserve the perfect (anti-)correlations for 
= and tt. The following lemma gives a sufficient condi- 
tion under which such a mixture can indeed give the desired 
cosine correlation 112 ill. 



Lemma 1. Let E((f>) be a (2tt -periodic, tt -anti-periodic, even) 
real function with a Fourier decomposition of the form 



(3) 



n>Q 



such that 



e\ > 0, e2n+i _ for all n > 1, 

£(0) = En>O e 2n+l = l, (4) 

^"(o) = -E„>o( 2 " + 1 ) 2e2 «+i <°- 

Then cos can be decomposed as 

cos = ^2 P2 m +iE((2m + 1)0) , (5) 

m>0 
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with P2m+i > for all m > 0. 

In particular, for = 0, one gets J2 m >oP 2m + 1 = ^ The 
coefficients P2m+i > can thus be interpreted as probabil- 
ities, and the kind of "inverse Fourier decomposition" (|5) is 
indeed a probabilistic mixture of correlations E((2m + 1)0) . 

A proof of Lemma 1 is given in Appendix B together with 
the explicit form of the P2m+i- It is easy to check that E\ 
satisfies the conditions Hence, there exist coefficients 
P2m+i > such that 2 m>0 P2m+i = 1 and 



P2r 



L .Ei((2m + 1)0) = cos0. 



(6) 



m>0 



Consequently, and since (2m + 1)0a + (2m + 1)4>b + (2m + 
l)4>c = (2m + 1)0, the following protocol gives the desired 
cosine correlation and solves our problem, with the same 3 
classical bits as in Protocol 1 : 

Protocol 2. Let Alice, Bob and Charlie share, in addition to 
the randomness already introduced in Protocol 1, a random 
variable M that takes the value M = 2m + 1 with proba- 
bility P2m+i, where {p2m+i}m>o ore the coefficients of the 
decomposition (|6]). 

After reception of their measurement settings 4>a,4>b and 
4>c, the three parties run Protocol 1 with input angles (2m + 
1)4>A, (2m + 1)0b and (2m + l)0c> respectively. 

Variants of our communication protocol. For conve- 
nience, let us from now on consider the equivalent 0/1 bit val- 
ues corresponding to the 3 parties outputs in Protocol 1 (or 2): 
a = , b = and c = ; the additions below will 



1=2 }j _ 
2 ' 2 2 

be modulo 2. Writing explicitly function 
of the classical communication (the bits tq, r b , t c ) that Alice 
receives, one has 



a T ii — a T oo + 1 an d a T0 oi — a 



TO 10 



(7) 



One can see that our communication protocol can actually 
be declined in different forms. In particular, Alice might not 
need to know the individual values of the bits r b and r c , but 
only their sum r bc = r b + t c . Charlie's bit t c could for ex- 
ample be sent to Bob instead; Bob would then send r bc to 
Alice, who would output a' TaTbc = a TO T bc o\ m the case when 
r b = t c = 1, the '+1' term in (0 can be introduced by Bob 
instead, who should output b' T = b + r b T c . This thus induces 
a protocol 1 ' , summarized as 



Prot. 1' 



C 

c' = 



A 



c, b> 



TbT c 



a T T ba ■ 



With similar considerations, one can come up with many 
different variants with varied communication patterns, such 
as, for instance: 



Prot. 1": 



Prot. V 



C 



B, B 



A^B 



C"=C, 0," = O to0 0, b"= b + T a T b +T a T c + T b T c 



B^C ^ A 



B 



c+T b T c , a'" bc = ao Tbc o, K' a = b+T T a , 



with r a = a To00 + a ToW and r a = a 0rbc o + a lTbc o- These 
variants and the original protocol look different, though they 
all require 3 bits of communication and lead to the same cor- 
relation. All of them have severe timing constraints (which 
is common for communication protocols): there are always 
some players that can't produce their output before some other 
partners receive their input and send them some information. 

Simulation with PR boxes. An interesting alternative to 
measure nonlocality is to estimate the number of nonlocal PR- 
boxes l22ll (some kind of "unit of nonlocality" Il23l0 required 
to simulate the correlations. Since the correlations we con- 
sider in this letter have no single- nor bi-partite marginals, all 
the variant communication protocols introduced above can be 
translated into PR-box based protocols 02 311 . Indeed, using (O 
for the original version of Protocol 1, one can always decom- 
pose the sum a + b + c as follows: 

ar a T b T c +b + c = a oao + b + c + T a (a aQa + a wa ) 

+ T b (aooo + aoio) + T T b (aooo+aoio + aioo + aiio) + UT C 
+ T c (aooo+aoio) + Tor c (aooo+aoio+aioo+aiio) ■ (8) 

The product terms in dS) can be generated by using nonlo- 
cal boxes: 5 PR-boxes can be used for the first 5 products (3 
between Alice and Bob, 1 between Bob and Charlie and 1 be- 
tween Alice and Charlie); the last product can be generated by 
a 3 -partite GHZ-box, which can in turn be constructed from 3 
PR-boxes 12411 . Hence, a total of 8 PR-boxes suffices to simu- 
late the tripartite GHZ correlations. 

Interestingly, the variant communication protocols de- 
scribed above all lead to a PR-box based protocol with the 
same configuration of 8 PR-boxes, all used precisely in the 
same way. In addition to this invariance, and similarly to 
quantum correlations, the PR-box based protocol does not suf- 
fer from any timing constraint. Hence, it might be a more 
faithful tool to measure quantum nonlocality (at least, for cor- 
relations with vanishing marginals) - this question is quite 
general and would require further scrutinies beyond the scope 
of this letter. Note finally that reciprocally, simulating the PR 
boxes by communication gives a systematic way to generate 
different variants of our initial protocol, depending on which 
way we the communication goes. 

Detection loophole. Another interesting connection is be- 
tween our communication protocol and simulation models 
based on the detection loophole 11121 12511 . For this connec- 
tion let us start for instance from the last variant of the com- 
munication protocols. In the detection-loophole-based proto- 
col, T b , r bc and r a are 3 additional shared random variables 
and each player outputs a bit if and only if the appropriate r 
agrees with the bit he should send in the communication pro- 
tocol. Hence, using variant 1"' of our protocol, the detection- 
loophole-based protocol simulates the GHZ correlations with 
"detection efficiencies" of 50% for Alice, Bob and Charlie. 
Other variant protocols can lead to detection-loophole-based 
protocols with asymmetric detection efficiencies. 

Conclusion. We have proven that 3 bits of communication 
(or 8 PR-boxes) suffice to simulate 3-qubit GHZ equatorial 
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correlations; hence the nonlocality of these correlations is at 
most of 3 bits (8 PR-boxes). In the course of our derivation, 
we introduced a strategy to obtain a cosine correlation as a 
mixture of other ("harmonic") correlations, via Lemma 1, that 
we believe could be used in other contexts as well. 

In this letter we considered correlations with vanishing 
single- and bi-partite marginals. If one considers also mea- 



surements on the GHZ state out of the equatorial plane II26I1 . 
or if one considers other states such as biased GHZ-like states 
for instance, then the marginals will no longer be random, and 
simulating the entire probability distribution is likely to be sig- 
nificantly harder 1U5I1 . 

Two other important open problems are the questions of the 
optimality of our protocol and of its generalization to more 
parties. For 3 parties, since the GHZ correlations are truly 3- 
partite [27], a minimum of 2 bits is necessary to connect the 
3 parties. We could find a 2-bit protocol (Protocol 1, with- 
out step 0, see 1I20I1 ) that gives stronger correlations than cos <p 
and that can approximate it to a very good accuracy, but not 
perfectly. For the iV-partite case, it is easy to generalize pro- 
tocol 1, again without step 0, using (N — 1) log 2 (iV — 1) bits 
of communication: divide the equator of the Bloch sphere into 
2(7V— 1) equal sectors, let each of the N— 1 last parties share a 
random angle tpi with Alice, and tell her in which sector their 
angle 4>j = cf)j + tpi (modulo ir) is. This leads again to a proto- 
col giving stronger correlations than cos (f> (actually, stronger 
and stronger as N increases), with a number of bits that is 
asymptotically equivalent to the lower bound derived in II 1611 
for the simulation of GHZ correlations (N log 2 N — 27V). Un- 
fortunately, we did not find a generalization that would give a 
correlation satisfying the assumptions of Lemma 1, so that 
the exact cosine correlation could then be obtained as in Pro- 
tocol 2. 

These observations lead us to formulate the following ques- 
tion: should we understand a "stronger" correlation as be- 
ing "more non-local"? If our goal is to quantify the power 
of nonlocality as a resource for achieving some information 
processing task, then the next question follows: is there any 
(useful) task, for which a stronger correlation might actually 
be less powerful than a weaker one? If this is not the case, 
then one could be happy with simulation protocols that give 
stronger correlations than the desired ones, and for this opera- 
tional interpretation of the nonlocality measure, we could con- 
clude that the nonlocality of the 3-partite GHZ correlations is 
at most 2 bits (or 3 PR-boxes), and that of the TV— partite GHZ 
correlations is at most (TV — 1) \og 2 (N — 1) bits. 

Nonlocal correlations are fascinating. First, because they 
can't be simulated by mere shared local variables; next, be- 
cause even if finite communication is allowed, their simula- 
tion remains tedious and quite artificial. Hence, simulating 
in particular quantum nonlocal correlations with classical re- 
sources, like shared local variables and communication, looks 
in general extremely difficult. This underlines the power of 
nonlocal correlations. Yet, such simulations seem to give a 
good measure of nonlocality (whether we are interested in the 
exact simulation or in the "operational nonlocality" measure), 



possibly the best together with PR-box based simulations, and 
provide the only story that takes place in space and time about 
how they could occur. 
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Appendix A: Calculation of E\((f>) 

In this Appendix we analyze the correlation obtained with 
Protocol 1. We first introduce two useful lemmas. 



Lemma 2. Let a be a vector on the equator of S 2 , with 
azimuthal angle (ft a- Consider two random vectors Ai,A2 
uniformly distributed on the half sphere S + (a) = {A G 
S 2 I a ■ A > 0}, and define X = ~ 1+ ^ 2 . Then the azimuthal 

angle p G [(ft a — f i 4>A + §] of X is distributed according to 



p{v) = o cos (^ ~ M- 



(9) 



Alternatively, if one doesn't want to restrict a priori ip to be 
in [4>a — § , (ft a + §], one can write, for any ip G [0, 2ir] (or 
any real interval with amplitude 2if), 



p(cp) = max(0, - cos(ip - <ft A ))- 



(10) 



Proof of Lemma [2} 

Denote by p$2 and p s + the uniform distributions on S 2 
and S + (a), respectively. For another vector b on the equator 
of S 2 , with azimuthal angle <ft B , let's calculate: 



Now, the term sign (b- (Ai +A2)) in the definition of I above 
is simply equal to sigii(cos(<^ — 0b)), so that the integral is 
also equal to 



r<t>A+- 

1 = cos{(ft B -(ftA)= I p{p) dp sign(cos(< ( 9 - <ft B )) ■ 

(ID 

Defining R(<ft) = ft* p{<p) dp (with R(<ft A + f ) = 1 due 
to normalization), the explicit calculation of ( fTTT i leads to 

7T 

cos(<ft B -4>A)= l-2R((ft B - -) if 4>A < (fts < <f>A + K, 

IT 

COs((ft B -(ftA)= 2R((ft B + if (ft A ~TT < 4>B < <ft>A ■ 

With now p = <ft B =F ^, both equalities above give 

7T 7T 

2R(<p) - 1 = sm(p - (ft A ) for (ft A - -z < P < 4>a + tt • 

After differentiation, one obtains 
1 



□ 



Lemma 3. After step of Protocol 1, pt G [0,27r] /* dis- 
tributed according to 



P(<Pb) = 7|sin2(^ 6 + <f> B )\ 



(12) 



p 5+ (Ai)dAi / / 9 i 5 + (A 2 )dA 2 sign(6-(Ai + A 2 )) 

S+(a) JS+(a) 

= Ps 2 i^i)dX 1 p S 2(X 2 )dX 2 sign(6-((TiA 1 + cr 2 A2)) 



p 5 2(Ai)dAi / /o 5 2(A 2 )dA 2 a x sign(6-(Ai + cricr 2 A 2 )) 
s 2 Js 2 



with &i = sign(a • A^). The last integrals were calculated 
by Toner and Bacon in lfl3ll . and were found to be equal to 

a ■ b = cos((ft B — 4>a)- 



Proof of Lemma \3} 

We now use the notation cr^ = sign(& • Aj). 

Suppose first that o\ = +1. Then Ao = Ai + cri(T 2 A 2 = 
o\X\ + cr 2 A 2 is the sum of two vectors uniformly distributed 
on S + (b). From Lemma [2] the distribution of its azimuthal 
angle (p G [0, 27r] is 

1 7T 

P+(fo) = max(0,-cos(<^ - -20b))) 
= max (0, - sm(p + 2<ft B )) . 



6 



In the case where o\ = — 1, Ao = Xi + 0102X2 = — ciM — 
(T2A2 is the sum of two vectors uniformly distributed on 
S + (— b). The distribution of its azimuthal angle tpo G [0, 2ir] 
is 

P-(<fo) = max (0, - sm(tp + 24>b + n)) . 

As Alice and Bob ignore the individual value of u\ (= ±1 
with equal probabilities), the overall distribution of ipg is 

Po(fo) = ^(p+{<po) +P-{<Po)) = j\sin((f Q + 2(j) B )\. 
The distribution of ip° = £ [0, n] is then 



It is convenient to use the Fourier decomposition 

sign(sina;) = | E«>o 5^h sin (( 2n + t0 calculate 
the integrals. One then easily gets 



n>0 V 



1 



(2n+l) 2 4-(2n+l) 2 



cos 



((2n+l)0) 



(13) 



One can finally check that Ei(<p) can a l so be written as 
£i(0) = 1 - ^"f 12 ^ for e [0,tt], from which the full 
function can be obtained by symmetry and periodicity. 



plivt) = 2po(2^) = 2 I MM + 20 s )| , 

and after adding £71-, with £ e {0,1} random, the distribution 

of ip b = ipl + £7r e [0, 27r] is finally 



P(<Pb) = 1 1 sin(2^ h + 20 B ) I ■ 

□ 



Let us now calculate the correlation obtained with Proto- 
col 1 . One can, for simplicity, directly integrate over the vari- 
ables 4>b = 4>b + fb an d 4>c = 4>c + fc\ from Lemma [3] 
4>b is distributed according to p(4>b) = j\ svh(2<Pb)\, while 
0c is uniformly distributed on [0, 27r]. One can easily check 
that the single- and bi-partite marginals vanish; the tripartite 
correlation writes 



{ah) 



2 71 



271 



0b 7 1 sin(20 B ) / sign(sin0 B ) sign(sin C ) 

4 Jo 2lr 

x sign (sm(4> B + 4>c - - n (4>b ) 77 - T c (0c ) ^ ) ) , 



Appendix B: Proof of Lemma 1 



Assume that E((j>) satisfies the assumptions (QUi. Define 



pi = — and, for m > 1, 



P2m+1 



_^ m — 1 m 

— E 52p 2fc + ie2n + i( 

ei ^— ' 



fc=0 n=l 



(2fe+l)(2n+l),2m+l 

(14) 



(with <5; ,j = 1 if i = j, Si.j = otherwise). 

Note that these coefficients P2m+i are non-negative. We 
will show that they lead to the decomposition (0). The proof 
is partly inspired by that of Lemma 3.1 in lE9ll : we divide it 
into 3 steps. 



• Step 1: We first prove that the coefficients P2m+i can 
be written as 

P2m+i -jLw.,... t»w... l~srj l"~J (15) 



and only depends on = 0a + s + Re- 
using the periodicity of the integrand function, one obtains 



where the (finite) sum is taken over all non-negative integers 
£3, £5, . . . such that 3' 3 x5 l5 x--. = 2m + 1. 



£x(0) := (a/3 7 ) = - 



6 C sin(20 B ) 



x sign( sin(0 s + 0c - 0)) ■ 



We prove it by induction. Note that (TT3T > holds for m = 0; 
suppose P2fe+i can be written as in (TT~5T > for all k < m. Then, 
with the notation i|j meaning "i divides j", one can write 



P2m+1 



1 

ei 



E 



e2n+l P^i±l 

2n + l 



s.t. 2n+l|2m+l 
1 m / 

d ^ I ei 



71=1, 

s.t. 2n+l|2m+l 



E 



(4 + 4 + ■■■)! 
• • • 



£3 
ei 



£5 
ei 



s.t. 3 e 35 f 5 



1„eIe _2m+l 



7 



P2m+1 — 



ei 
1 

ei 



E 



E 



(4 + + ...)' f e 3 



n=i, «3.-*Sn+l.- 
s.t. 2n+l|2m+l g t 3 * 3 . . . (2n+l) f 2r> + l + i . . . =2m+ i 



E E 

71=1, « 3 ,f 5 ,... 

s.t. 2n+l|2m+l s .t. 3'35 i 5 — =2m+l 



^2n+l 



(4 + 4 + •■•)! f e 3 



ei 



e2ra+l 

ei 



i+i 



t3 T t 5 



+ . . . £ 3 !4! • 



£5 

ei 



where we relabeled (^2«+i + 1) — > £271+1- After exchanging the two sums, and using the fact that ErILi e 3 +e^+ ~ ^' one 
obtains (l~5[ as desired. 



• Step 2: We now show that E-P2m+i converges abso- 
lutely. 

Eq. $15[ can be written as 



P2m+1 — 



+ 1 e x (2m+l) 



l i 

4F ^£ 3 ,£ 5t ... 



3+^5 + -)! 



3/2 \ f 3/ 3/2 

ei 



All the terms in the sum are non-negative. If one extends the 
sum to all integers 4>4j • • • > one gets the upper bound 



P2m+1 < 7" 



L! 



ei (2m+l) 3 / 2 Z^L=0 l^l z ,l b ,... l s \t 5 \... 



where in the second sum, the indices are now such that 4 + 
4 + • • ' = L. 

From the assumptions ©, E(2™ + 1) 3 ^ 2 (— e 2 „+i) < 
E(2« + 1) 2 (— e 2n+ i) converges absolutely, and one can ap- 
ply the multinomial theorem to calculate the inner sum; one 
finds that this sum is (C3/ 2 ) L , where C3/2 = EriLi ( — 

(2n+l) 3/2 e 2 „ + i \ 

Now'o < E,T=i(2^ + l) 3/2 (-e 2 „ +1 ) < £~ =1 (2n + 
1) 2 (~ e 2n+i) < ei by assumption 1 ; hence < C 3 / 2 < 1, and 
therefore one gets 

< P2m+1 < 57 ( 2m +l)3/2 ElU ( C '3/2) L 
_ 1 1 1 



ei l-C 3/2 (2m+l) 3 / 2 ' 

which implies that EP2m+i converges absolutely. 

Note also that the assumptions (01 imply that E e 2n+i con- 
verges absolutely as well. 



• Step 3: Conclusion 

The fact that Ef>2m+i an d E e 2«+i are absolutely conver- 
gent allows one to calculate the following infinite sums: 



P2m+iE((2m + l)<j>) 



m>0 



m>0 n>0 
fc fe 

E(EEp 2m + ie2) 

fe>0 m=0n=0 



+1 <5(2m+l ) (2«+l ) , 2fc+l) C0S ( ( 2 fc + 1 ) </>) 



By definition (TBI) , the double sum inside the brackets is equal 
to Sk,i, and one obtains, as desired, 

P2m+iE((2m + 1)0) = cos(0) . (16) 



m>0 



□ 



The strict inequality in what precedes should actually be replaced by an 
equality in the case where ei = 1 and e2n+i = for all n > 1. But the 
conclusion still holds in that trivial case. 



